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I.  Introduction  and  Program  Objective 

Current  methods  of  lifetime  assessment  leave  much  to  be  desired.  Typically,  the  expected  life 
of  a  full-scale  component  exposed  to  a  complex  environment  is  based  upon  empirical 
interpretations  of  measurements  performed  on  microscopic  samples  in  controlled  laboratory 
conditions.  Extrapolation  to  the  service  component  is  accomplished  by  scaling  laws  which,  if 
used  at  all,  are  empirical;  little  or  no  attention  is  paid  to  synergistic  interactions  between  the 
different  components  of  the  real  environment.  With  the  increasingly  hostile  conditions  which 
must  be  faced  in  modem  aerospace  applications,  improvement  in  lifetime  estimation  is  mandated 
by  both  cost  and  safety  considerations. 

This  program  aims  at  improving  current  methods  of  lifetime  assessment  by  building  in  the 
characteristics  of  the  micro-mechanisms  known  to  be  responsible  for  damage  and  failure.  The 
broad  approach  entails  the  integration  and,  where  necessary,  augmentation  of  the  micro-scale 
research  results  currently  available  in  the  literature  into  a  macro-scale  model  with  predictive 
capability. 

In  more  detail,  the  program  will  develop  a  set  of  hierarchically  structured  models  at  different 
length  scales,  from  atomic  to  macroscopic,  at  each  level  taking  as  parametric  input  the  results 
of  the  model  at  the  next  smaller  scale.  In  this  way  the  known  microscopic  properties  can  be 
transported  by  systematic  procedures  to  the  unknown  macro-scale  region.  It  may  not  be  possible 
to  eliminate  empiricism  completely,  because  some  of  the  quantities  involv^  cannot  yet  be 
estimated  to  the  required  degree  of  precision.  In  this  case  the  aim  will  be  at  least  to  eliminate 
functional  empiricism.  Restriction  of  empiricism  to  the  choice  of  parameters  to  be  input  to 
known  functional  forms  permits  some  confidence  in  extrapolation  procedures  and  has  the 
advantage  that  the  models  can  readily  be  updated  as  better  estimates  of  the  parameters  become 
available. 


II.  Program  Organization 

The  program  has  been  organized  into  specific  tasks  and  subtasks  as  follows. 

Task  100.  Lifetimes  of  metallic  dispersed-phase  composites 

Most  service  materials  fall  into  the  category  of  dispersion-hardened  metallic  composites.  This 
task  will  consider  the  problem  of  dispersion  hardened  materials  in  general,  but  with  two  specific 
materials,  NiAl  and  MoSi2/SiC  in  mind. 

Task  110.  Identincation  and  modelling  of  micromechanisms 

The  purpose  of  this  task  is  to  determine  what  micromechanisms  are  operative  in  the  high- 
temperature  deformation  of  dispersion-hardened  materials.  In  the  general  case  this  will  be  done 
by  a  literature  search.  For  specific  materials,  the  micromechanisms  will  be  determined  from  the 
experimental  program  at  NRL.  Once  identified,  each  of  these  micromechanisms  will  be 
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modelled,  in  order  to  determine  what  are  the  critical  parameters  which  determine  its  effect  on 
plastic  flow  and  values  for  these  parameters.  Also  to  be  determined  is  whether  the  modelled 
critical  values  are  dependent  on  quantities  which  must  be  obtained  from  a  smaller  scale  model. 

Task  111.  Equiaxed  dispersoids 

This  task  will  consider  dispersions  of  the  type  encountered  in  NiAl-like  materials.  That  is,  the 
dispersoids  are  considered  to  be  small  compared  to  the  grain  size.  The  term  ’equiaxed’  is  used 
because  the  particles  are  roughly  of  the  same  size  in  all  three  dimensions.  However,  this  is  not 
a  requirement  for  this  task.  Rather,  it  is  necessary  that  the  particles  not  be  too  large  in  the 
dimension  normal  to  the  slip  plane,  so  that  they  can  be  surmounted  with  relative  ease  by  cross¬ 
slip  and/or  climb  without  the  generation  of  appreciable  back-stress. 

Task  112.  Anisotropic  dispersoids 

This  task  covers  the  case  of  dispersoids  which  are  elongated  in  the  direction  normal  to  the  slip 
plane.  An  example  is  SiC  fibers  in  MoSi2.  In  this  case,  plastic  flow  around  the  dispersoids  takes 
place  by  a  combination  of  glide  and  climb,  but  is  a  protracted  process  during  which  large 
stresses  acting  in  opposition  to  the  applied  load  are  developed. 

Task  113.  Grain  boundary  effects 

This  task  will  examine  the  role  of  grain  boundary  processes  in  high-temperature  deformation. 

Task  120.  Macroscopic  stochastic  model  for  creep 

In  real  materials  it  is  likely  that  more  than  one  mechanism  will  be  operative,  either  in  parallel 
or  in  series.  The  information  gained  in  task  1 10  is  not  sufficient  to  describe  this  situation.  Once 
the  critical  parameters  for  individual  mechanisms  have  been  determined,  it  is  necessary  to 
combine  them  in  a  macroscale  stochastic  model.  This  will  be  done  by  determining  critical 
stresses  and  activation  enthalpies  as  a  function  of  local  geometry  and  using  these  values  in  a 
finite-temperature  simulation  of  creep  through  a  random  array  of  dispersoids.  Careful  attention 
must  be  paid  to  possible  interactions  between  mechanisms. 

Task  130.  Extension  to  cyclic  deformation 

The  final  step  in  task  100  is  to  extend  the  results  to  the  case  of  cyclic  deformation.  Irreversibility 
is  an  intrinsic  feature  of  the  model  in  task  120.  However,  it  is  likely  that  other,  as  yet 
unrecognized,  characteristics  of  cycled  deformation  will  have  to  be  considered. 

Task  200.  Lifetimes  of  piezoelectric  ferroeiectrics 

Failure  in  cyclic  loading  of  sensors  and  actuators  formed  from  lead  zirconate  titanate  (PZT)  is 
a  continuing  problem.  PZT  is  a  ceramic  and  therefore  differs  from  the  materials  considered  in 


3 


task  100  in  that  plastic  deformation  is  not  involved.  This  task  will  examine,  modelling  as 
necessary,  the  operation  of  PZT  devices,  in  order  to  determine  the  factors  governing  lifetime 
limitation. 

Task  300.  Reporting 

Running  concurrently  with  tasks  100  and  200,  this  task  will  inform  the  Navy  Program  Manager 
and  Contracting  Officer  of  the  technical  and  fiscal  status  of  the  program  through  R&D  status 
reports. 
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III.  Technical  Progress 
Task  120.  Macioscopic  stochastic  model  for  creep 

A  numerical  model  has  been  developed  to  predict  creep  strengths  and  lifetimes  for  a  metallic 
matrix  with  an  arbitrary  dispersion  of  seco^-phase  particles  of  varying  size  and  strength.  An 
analytic  treatment  of  the  limiting  case  of  impenetrable  particles  is  given  in  the  preprint 
reproduced  in  Appendix  A. 

a.  Theory 


a.l  Description  of  Model 

Figure  1  defines  the  parameters  which 
describe  the  geometry  of  the  general 
dispersion.  Dispersoids  (assumed  to  be 
spherical)  intersect  the  slip  plane  in 
circles  (shown  in  red  in  Figure  1)  of 
different  diameters,  Dj,  and  are 
separated  by  distances  Ljj.  The  angle 
between  the  dislocation  arms  at  the  ith 
dispersoid  is  6^.  The  value  assumed  by 
the  angle  6  at  the  point  of  Orowan  looping  defines  the  strength  of  the  dispersoid,  as  described 
in  the  previous  report'.  An  important  feature  of  the  model  is  that  the  mechanical  properties  of  the 
dispersion,  both  athermal  and  thermal,  are  defined  completely  by  the  parameters  L,  D  and  0'.  To 
recapitulate,  L  and  D  define  the  line  tension,  T,  of  the  dislocation,  via  the  relations 

T  =  A-^lnJT  1=1+1  (1) 

4x  X  L  D 

where  p  is  the  shear  modulus,  b  is  the  magnitude  of  the  Burgers'  vector  and  A  is  a  constant  with 
the  value  1  for  edge  dislocations  and  1/(1 -v)  for  screw  dislocations  (v  is  Poisson's  ratio).  The 
strength  of  the  dispersoid  is  determined  by  the  angle  0,  by  way  of  the  expression 

xhL  =  22l;os-|  (2) 


Figure  1  Definition  of  the  geometric  disposion 
parameters  D,  L  and  6 


in  which  t  is  the  effective  stress.  The  left  hand  ade  of  (2)  represents  the  forward  force  on  the 
dislocation,  while  the  right  hand  side  is  the  opposing  force  exerted  by  the  dispersoid. 

The  formalism  of  (1)  and  (2)  follows  the  generalization  by  Duesbery  and  Sadananda^  of  a  model 
put  forward  by  Bacon,  Kocks  and  Scattergood^  and  has  b^n  shown  to  be  compatible  with  a  wide 
range  of  experimental  observations'. 
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Key  innovations  of  the  numerical  model  are  as  follows. 

★Explicit  temperature  dependence  in  both  strength  and  creep  calculations  is 
incorporated.  The  activation  enthalpy  for  any  dispersoid  can  be  calculate  by 
integrating  the  expression  (2)  with  respect  to  6  from  the  equilibrium  to  the  critical 
value. 

★The  estimation  of  lifetimes  of  order  up  to  years  from  a  model  with  a  time  scale 
of  picoseconds  is  impractical.  A  novel  algorithm  has  been  developed  which 
permits  prediction  of  long-term  creep  properties  from  a  microsc^e  model. 
Knowledge  of  the  activation  enthalpy  for  each  dispersoid  along  the  dislocation 
line  permits  computation  of  the  probability  of  activation.  Combining  this  with  the 
vibration  frequency  of  the  dislocation  allows  determination  of  the  mean  waiting 
time  before  activation.  The  algorithm  determines  the  dispersoid  with  the  least 
waiting  time.  This  dispersoid  is  then  activated,  the  elapsed  time  is  increased  by 
an  increment  equal  to  this  waiting  time  and  the  wait  time  for  all  other  dispersoids 
is  decreased  by  the  same  increment.  Any  true  instabilities  are  allowed  to  take 
precedence  over  thermal  events. 

★Arbitrary  combinations  of  obstacle  strength  and  size  are  permitted.  The 
intersection  of  each  dispersoid  with  the  slip  plane  is  chosen  to  lie  at  a  random 
position  along  the  dimension  of  the  dispersoid  normal  to  the  slip  plane. 

★Intersection,  Orowan  looping  and  diffusive  climb  mechanisms,  together  with 
thermally  activated  processes,  are  all  permitted  and  form  competing  mechanisms 
in  the  model.  Climb  is  simulated  by  assigning  a  wait  period  equal  to  the  time 
required  for  Einsteinian  diffusive  drift  over  the  shortest  distance,  normal  to  the 
slip  plane,  to  the  edge  of  the  dispersoid. 

★Duesbery-Sadananda  (DS)  line  tension  formalism  is  incorporated  in  the 
harmonic  mean  approximation. 


a.2  Strength  of  Random  Dispersions 

The  model  can  be  used  to  calculate  the  athermal  strength  of  random  dispersions.  This  is  of 
interest  because  it  permits  comparison  of  the  predictions  of  the  model  with  those  of  earlier, 
simpler  approaches.  Figure  2  shows  the  computed  strengths  of  random  dispersions  as  a  function 
of  critical  breaking  angle  for  a  constant  number  density  of  dispersoids  (i.e.  constant  mean 
separation)  and  three  different  dispersoid  diameters.  The  unit  of  strength  is  the  strength  of  a 
square  array  of  point  obstacles  in  the  classical  Friedel  line  tension  approximation,  that  is,  /tb/L. 
This  choice  permits  the  primary  strength  dependence,  on  L  ',  to  be  factored  out,  thereby 
emphasiizing  the  influence  of  the  different  line  tension  models.  The  black  curve  in  Figure  2 
shows  the  behavior  predicted  by  the  current  computer  codes  for  Friedel  line  tension.  This  agrees 
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Figure  2  Strength  of  random  dispersions  for  different  dispersoid  sizes 


closely  with  the  early  calculations  of  Foreman  and  Makin^  and  serves  to  confirm  the  validity  of 
the  model.  The  remaining  curves  show  the  behavior  for  DS  line  tension  for  dispersoid  diameters 
of  100b  (blue  curve),  20b  (green)  and  4b  (red).  Common  to  all  the  DS  curves  is  a  strength  of 
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30%  or  more  greater  than  the  classical  value.  This  is  most  noticeable  for  large  breaking  angles 
(weak  dispersoids)  and  small  diameters.  In  practical  terms,  this  means  that  for  small  dispersoids, 
the  intrinsic  strength  of  the  dispersoids  is  not  of  great  consequence,  provided  that  the  breaking 
angle  is  less  than  about  120°.  Figure  3  shows  a  similar  strength  plot,  for  a  constant  dispersoid 
diameter  of  50b  and  for  mean  separations  of  810b  (red),  362b  (green)  and  162b  (blue).  Once 
again,  it  is  evident  that  the  more  accurate  DS  line  tension  predicts  strengths  considerably  in 
excess  of  that  of  the  classical  model. 

a.3  Creep  Lifetimes 

The  main  purpose  of  the  current  model  is  the  prediction  of  creep  lifetimes.  Using  the  novel 
algorithm  described  above,  this  has  been  done  for  a  representative  dispersion  with  volume 
fraction  of  5%  and  dispersoid  diameter  of  100b.  In  a  random  arrangement,  these  values 
correspond  to  a  mean  apparent  separation  <L>  =  324b  and  a  mean  intersected  diameter  <D>  = 


Figure  4  Creep  lifetime  as  a  function  of  aeep  stress  ratio  for  a  self-difliision  eno^gy  of  3eV 

(29iKJ/mol) 


78b.  A  breaking  angle  of  90°  has  been  chosen  as  typical  of  many  experimental  measurements  and 
theoretical  simulations'.  In  Figure  3  is  plotted  the  quantity  log,o(pb^A/ec)  (which  will  be  called 
the  lifetime  factor),  where  p  is  the  dislocation  density  and  A  is  the  lifetime  in  seconds  to  a 
critical  strain  of  e^,  against  the  ratio  of  the  creep  stress  to  the  athermal  strength.  The  particular 
values  in  Figure  3  correspond  to  a  ratio  of  p/e,.  =  10'®,  as  would  be  the  case,  for  example,  for 
a  dislocation  density  of  10*  cm  ^  and  a  critical  strain  of  0.01.  It  can  be  seen  from  Figure  3  that 
the  model  can  predict  lifetimes  over  a  range  of  20  orders  of  magnitude.  Each  point  is  an  average 
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of  five  calculations  based  on  different  random  number  sequences,  and  each  calculation,  taking 
15-30  minutes  on  a  fast  personal  computer,  samples  the  creep  motion  of  a  single  dislocation  over 
a  distance  of  lO’b  in  a  slip  plane  segment  of  width  lO'b.  The  three  curves  in  figure  3  refer  to 
temperatures  of  500“C  (red),  1000“C  (green)  and  1500“C  (blue).  There  are  three  clear  regimes. 
For  high  creep  stress  ratios,  of  order  0.75,  the  lifetime  is  very  short  For  creep  stress  ratios  0.65 
to  0.5,  the  lifetime  increases  very  rapidly,  finally  saturating  into  a  low-index  power  law 
dependence  for  creep  stress  ratios  of  less  than  0.45.  The  physical  reasons  for  the  three  stages  are 
as  follows.  At  high  creep  stress  ratios,  the  inter-obstacle  spacing  is  given  accurately  by  the  mean 
apparent  separation,  and  each  successful  thermal  activation  results  in  10^  -  10^  secondary, 
mechanical  activations  (i.e.  induced  instabilities).  In  the  transition  region,  the  number  of 
secondary  activations  drops  rapidly.  At  the  same  time,  the  dislocation  curvature  decreases; 
consequently  the  true  separation  of  the  obstacles  along  the  dislocation  line  increases,  so  that  there 
are  fewer  activation  sites.  In  the  low-stress  limit,  the  number  of  secondary  activations  saturates 
at  <  1  and  the  true  separation  of  obstacles  approaches  3<L>. 

These  factors  alone  are  insufficient  to  explain  the  curves  in  Figure  3  completely.  The  activation 
energies  for  overcoming  the  dispersoids  by  thermally  activated  glide  can  be  very  large,  especially 
at  low  stresses.  With  glide  processes  alone,  the  lifetimes  become  extremely  large.  Observation 
of  the  development  of  the  substructure  indicates  that  the  dislocations  become  ’stuck’  on 


Figure  5  Creep  lifetime  as  a  function  of  creep  stress  ratio  for  a  self-diffusion  energy  of  2eV 

(194KJ/inol) 


sequences  of  "hard"  obstacles.  To  keep  the  lifetimes  at  reasonable  levels,  it  is  necessary  to  allow 
for  the  overcoming  of  dispersoids  by  climb.  In  Hgure  3,  the  self-diffusion  energy,  which  governs 
the  climb  process,  is  set  at  3eV  (291KJ/mol).  Figures  4  and  5  show  the  change  in  behavior  when 
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self-duffusion  energy  is  set  to  2eV  ( 194KJ/mol)  and  leV  (97K.!/mol),  respectively.  The  low-stress 
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creep  stress  ratio  x  10 


Figure  6  Creep  lifetime  as  a  function  of  creep  stress  ratio  for  a  self-diffusion  energy  of  leV 

(97KJ/mol) 

lifetime  factor  at  500°  C  drops  by  13  orders  of  magnitude  as  the  self-diffusion  energy  drops  from 
3eV  to  leV,  exactly  as  expected  from  a  mechanism  controlled  by  diffusion. 

b.  Comparison  with  experiment 

The  results  reported  above  are  "hot  off  the  press"  and  as  yet  there  has  been  no  time  to  make  the 
comparison  with  experiment. 

c.  ImpUcaUons 

★That  low -stress  creep  is  dominated  by  dislocation  climb  and  hence  by  the  self-diffusion  energy. 

★Long-term  creep  can  be  modelled  by  computer  simulation,  without  the  need  for  empirical 
extrapolation. 

d.  Plans  (given  continued  funding) 

★The  results  of  the  creep  model  will  be  compared  both  with  experimental  results  and  with  the 
standard  engineering  formulae  used  to  extrapolate  long-term  creep  properties*. 
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Task  200.  Lifetimes  of  piezoelectric  ferroelectrics 

The  utility  of  PZT  as  a  ’smart  material’  is  severely  limited  by  its  premature  mechanical  failure 
when  subjected  to  cyclic  loading.  The  cause  of  this  failure  forms  the  critical  question.  The 
central  problem  is  that  an  irreversible  mechanism  is  necessary  to  explain  the  experimental 
observations,  while  the  overt  physical  properties  -  that  is,  piezoelectricity,  ferroelectricity  and 
ferroelasticity  -  are  all  symmetric  to  reversal.  A  theory  will  be  presented  which  introduces  an 
intrinsic  irreversibility  and,  at  the  same  time,  explains  all  salient  experimental  observations. 

The  key  experimenal  observations  are  as  follows. 

★In  the  poled  state,  ferroelectric  domains  are  elongated  strongly  in  the  direction  of  the  poling 
field  and  are  separated  by  walls  with  little  curvature. 

★Under  cyclic  mechanical  loading  at  70%  of  the  monotonic  fracture  stress  (i40MPa),  the 
lifetime  is  in  excess  of  10^  cycles.  At  stress  amplitudes  from  75%  to  95%  of  the  fracture  stress, 
the  lifetime  decreases  sharply  from  10*  to  10^  cycles)*.  The  effect  of  temperature  is  unknown. 

★Under  cyclic  (resonant)  electrical  loading,  the  lifetime  depends  on  temperature.  For 
deformation  above  80®  C,  failure  is  much  more  rapid  than  below  this  crtical  temperature. 
Predeformation  below  80°  C  prior  to  high-temperature  deformation  leads  to  even  more  rapid 
failure.  Damage  is  not  detectable  for  deformation  below  80°  C,  but  is  observed  as  intergranular 
cracking  at  higher  temperatures*. 

a.  Theory 

a.l  Domain  Wall  Structure 

There  are  two  basic  classes  of  domain  wall  in  PZT  and  other  ferroelectrics  belonging  to  the 
perovskite  family.  The  180°  wall  separates  domains  in  which  the  polarization  is  reversed,  but 
the  lattice  remains  commensurate.  On  the  other  hand,  the  polarization  (which  is  always  parallel 
to  the  tetragonal  c-axis)  changes  by  90°  on  passage  through  a  90°  domain  wall,  so  that  there  is 
an  incommensuracy  where  the  c-axis  on  one  side  of  the  wall  meets  one  of  the  a-axes  on  the 
other.  In  strictly  geometric  terms,  the  90°  wall  is  an  interface  between  two  crystals  with  different 
lattice  constants.  This  is  similar  to  the  problem  of  an  epitaxial  layer  on  an  incommensurate 
substrate  and  the  physical  consequences  are  well  known'.  For  thin  layers  and  small  misfit,  the 
incommensuracy  is  accommodated  by  coherency  strains.  As  the  thickness  of  the  layer  (or  in  the 
case  of  a  ferroelectric,  the  domain  width)  increases,  it  becomes  increasingly  advantageous  to 
concentrate  the  misfit  into  dislocations.  With  an  approximate  theory,  Matthews’  estimates  the 
maximum  misfit  £  that  can  be  accommodated  by  a  layer  of  thickness  h  without  instability  toward 
dislocation  formation  as 
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so  that,  for  example,  dislocations  would  be  preferred  in  a  layer  of  thickness  100b  if  the  misfit 
strain  is  greater  than  0.0017.  Therefore 
for  PZT,  for  which  the  misfit  is  1  -  3%, 
depending  on  the  composition, 
dislocated  domain  walls  are  certainly  to 
be  expected.  An  example  of  a  domain 
wall  dislocation  is  shown  schematically 
in  Figure  7.  In  Figure  7  the  wall  is 
shown  as  a  dashed  line  and  the  PZT 
unit  cell  is  outlined  at  the  top  left  of  the 
figure.  The  dislocations  will  be 
regularly  spaced  along  the  domain  wall 
and  will  have  Burgers  vectors  which  lie 
in  the  plane  of  the  wall.  Conservative 
motion  will  be  possible  only  in  the 
plane  of  the  wall;  the  mobility  of 
dislocations  in  the  wall  is  uncertain,  but 
will  probably  be  low  and  thermally 
activated,  as  is  usual  for  ceramic 
materials. 

There  are  two  sub-species  of  90“  domain  wall.  These  are  shown  schematically  in  Figure  8,  which 
shows  unit  cells  of  the  structure  on  either  side  of  the  domain  wall  (green).  The  direction  of 
polarization  in  the  reference  domain  cell  is  indicated  by  a  red-tipped  arrow,  that  in  the 
neighboring  domain  cell  by  a  yellow-tipped  arrow.  The  tetragonality  is  exaggerated  for 
convenience.  The  two  diagrams  in  Figure  8  show  that  the  lattice  in  the  neighboring  domain  can 
lie  in  two  different  orientations.  In  the  upper  diagram,  the  polarizing  vector  lies  normal  to  the 
domain  wall  and  there  is  misfit  in  only  one  of  the  two  cell  dimensions  in  the  wall  plane.  A  single 
set  of  dislocations  suffices  to  relieve  this  misfit  strain.  In  the  lower  diagram,  the  neighboring 
polarization  vector  lies  parallel  to  the  domain  wall  and  there  is  misfit  in  both  cell  dimensions  in 
the  wall  plane.  To  relievce  this  strain,  two  sets  of  dislocations,  perpendicular  to  one  another,  are 
required.  The  description  can  readily  be  generalized.  The  upper  diagram  in  Figure  8  describes 
what  will  be  called  a  type  I  wall,  which  has  one  normal  polarization  vector.  A  type  II  wall  has 
two  parallel  polarization  vectors.  In  a  first  approximation,  the  type  II  wall,  with  double  the 
density  of  dislocations,  will  have  twice  the  energy  of  a  type  I  wall.  It  can  therefore  be  expected 
to  occur  only  when  geometrically  necessary. 

The  domain  wall  structure  is  directly  pertinent  to  domain  wall  motion,  and  therefore  a  few 
preliminary  remarks  are  in  order.  The  movement  of  a  wall  normal  to  itself  requires  transport  of 
misfit  strain.  If  the  misfit  has  been  concentrated  into  dislocations,  dislocation  climb,  a  non¬ 
conservative,  thermally  activated  process  involving  the  absorbtion  or  emission  of  vacancies,  will 
be  necessary.  In  this  case  wall  motion  will  be  relatively  difficult,  sensitive  to  temperature  and 


Figure  7  A  domain  wall  dislocation  in  PZT 
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deformation  rate.  The  climb  process  will  be  particularly  difficult  in  PZT,  because  the  vacancies 
must  be  emitted  in  pairs  to  maintain 
charge  neutrality. 

On  the  other  hand,  if  the  misfit  is 
sufficiently  small  that  coherent  walls 
are  favored,  wall  motion  should  be 
rapid  and  insensitive  to  temperature  and 
rate.  There  is  an  intriguing  possibility 
that  both  types  of  motion  may  occur. 

The  misfit  in  PZT  decreases  with 
increasing  temperature,  so  that  a 
transition  from  dislocation-pinned  walls 
to  freely  mobile  coherent  walls  is  quite 
possible.  A  transition  of  this  type  could 
explain  the  temperature  dependence  of 
cyclic  failure  in  electrical  cycling  noted 
above.  More  detailed  calculations  and 
critical  experiments  are  required  to 
pursue  this  hypothesis. 

a.2  Domain  Wall  Motion 

Because  domain  walls  in  ferro-electrics 
are  stable  (or,  at  least  metastable)  lattice  defects,  they  do  not  couple  directly  to  external  fields, 
instead  responding  to  thermodynamic  forces  -  that  is,  to  free  energy  gradients.  In  general,  both 
electrostatic  and  mechanical  terms  must  be  taken  into  account.  A  simplified  mean-field  theory 
will  be  presented.  Although  crude,  this  theory  demonstrates  that  there  is  a  fundamental 
irreversibility  which  will  lead  to  failure  in  cyclic  loading. 

There  are  two  sources  of  internal  energy  -  m^hanical  and  electrostatic.  Consider  first  the 
mechanical  energy.  The  initial  state  is  polycrystalline  cubic  PZT;  a  single  grain  will  be  treated. 
It  is  assumed  that  after  transformation  to  the  tetragonal  structure  the  stresses  within  the  grain  are 
uniform.  Then  if  a  volume  fraction  f  has  its  c-axis  lying  parallel  to  the  z-axis  of  an  external  frame 
of  reference  which  coincides  with  the  axes  of  the  original  cubic  structure,  the  elastic  energy  is 
given  by 


Figure  8  The  two  types  of  90"  domain  wall 


(l+3f)  (4) 

where  p  is  the  shear  modulus  and  c,  a  are  the  tetragonal  lattice  constants.  In  the  transformed  state 
it  is  expected  that  the  orientations  of  the  tetragonal  cells  will  be  random,  so  that  f  =  1/3  and  U^, 
is  a  maximum.  For  c/a  -1  =  6  =  0.017  (for  composition  close  to  the  morphotropic  phase  boundary 
and  a  temperature  of  298K’)  and  p  =  40GPa,  U,,  is  about  4.42MJ/m\ 
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After  poling  with  an  electric  Held  along  the  z  axis,  the  electrostatic  energy  can  be  estimated  from 
the  remanent  polarization,  P,,  which  has  a  value*  of  0.40C/m^  for  the  composition  and 
temperature  given  above.  Assuming  that  P,  is  proportional  to  5,  the  electrostatic  energy  is  given 
by 


= 


2  KK^ 


(5) 


in  which  k  is  the  dielectric  constant  along  the  z  axis  (  ••  1000)  and  is  the  permittivity  of  free 
space  (  =  8.8xl0  ‘^F/m).  The  maximum  electrostatic  energy,  from  (5),  of  9.0MJ/m^  is  comparable 
with  the  maximum  elastic  energy. 


The  two  terms  (4)  and  (5)  form  contributions  to  the  internal  energy  which  change  as  domain 
walls  move  and  therefore  exert  generalized  forces  on  the  walls,  as 


Pel 


df  dx 


dU„  df 
df  dx 


(6) 


where  x  is  measured  normal  to  the  domain  walls  and  the  poling  axis,  and  is  also  the  direction 
of  tensile  and  compressive  stress  during  mechanical  load  cycling. 

The  final  force  on  the  domain  walls  occurs  when  an  external  stress  is  applied.  This  is  seen  best 
by  reference  to  Figure  9,  which  shows  a  type  I  domain  wall  (vertical  line)  separating  domains 
with  cell  shape  and  polarization  shown  by  rectangles  and  red-tipped  arrows.  For  a  tensile  stress 
a  acting  normal  to  the  domain  wail,  a  movement  of  the  wall  to  the  left  through  a  distance  dx 
causes  the  external  stress  to  do  work  -o5dx.  Therefore  there  is  a  configurational  force  of  -o8  on 
the  wall,  as  shown.  Figure  9  also  shows  the  sense  of  the 
forces  Fj,  and  F„.  For  a  compressive  applied  stress  the  force 
aS  is  reversed  in  sign,  but  the  other  forces  remain 
unchanged. 

The  essential  irreversibility  of  domain  wall  motion  stems 
from  the  imbalance  of  F^,  and  F„.  The  resultant  of  these 
two  forces  maintains  a  constant  bias  while  the  external 
stress  reverses  during  cyclic  deformation.  Converting  all 
forces  (per  unit  area)  to  units  of  fi,  for  ease  of  comparison 

—  =  1.506^(f--^)  -  1.566*r  T  6-^  (7) 

|i  3  |i 

in  which  the  terms  on  the  right-hand  side,  in  order,  are  due 
to  the  elastic  energy,  the  electrostatic  energy  and  the 
applied  stress.  The  sign  of  the  latter  term  is  to  he  taken  as 
negative  for  tension,  positive  for  compression.  Effectively 
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Figure  9  The  forces  acting  on  a 
domain  wall 
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separated  by  walls  with  little  curvature. 

A  domain  wall  which  is  normal  to  the  electric  field  has  a  surface  charge  which  experiences  a 
force  due  to  the  field  which  promotes  motion  of  the  wall  in  the  direction  of  the  field  axis.  Walls 
interacts  directly  parallel  to  the  field  experience  no  such  force.  Hence  domains  tend  to  elongate 
parallel  to  the  field  axis,  growing  until  terminated  by  the  crystal  surface  or  a  grain  boundary. 
The  remaining  walls  are  constrained  by  misfit  dislocations  to  maintain  long-range  linear  form. 
At  the  short  range,  curvature  of  the  walls,  which  will  be  low-energy  type  I  walls,  requires 
formation  of  small  segments  of  type  n  walls,  and  is  therefore  inhibited. 

MJrtder  cyclic  mechanical  loading  at  70%  of  the  monotonic  fracture  stress  (140MPa),  the 
lifetime  is  in  excess  of  l(f  cycles.  At  stress  amplitudes  from  75%  to  95%  of  the  fracture  stress, 
the  lifetime  decreases  sharply  from  l(f  to  l(f  cycles  f.  The  effect  of  temperature  is  unknown. 

As  the  crystal  depoles  by  domain  wall  movement,  the  domains  lengthen  in  a  direction  parallel 
to  the  stress  axis  and  shorten  in  the  poling  direction.  This  creates  cracks  at  the  grain  boundaries 
which  lengthen  systematically  during  cyclic  stressing,  due  to  the  intrinsic  irreversibility  of 
domain  wall  motion,  until  a  crack  of  critical  length  is  created  and  fracture  occurs.  This  process 
is  strongly  dependent  on  temperature  and  stress,  because  the  controlling  mechanism  is  dislocation 
climb. 

*Under  cyclic  (resonant)  electrical  loading,  the  lifetime  depends  on  temperature.  For 
deformation  above  Sff’  C,  failure  is  much  more  rapid  than  below  this  critical  temperature. 
Predeformation  below  8(P  C  prior  to  high-temperature  deformation  leads  to  even  more  rapid 
failure.  Damage  is  not  detectMe  for  deformation  below  8&  C,  but  is  observed  as  intergranular 
cracking  at  higher  temperature 

There  is  a  critical  temperature  above  which  dislocations  will  be  expelled  from  the  domain  walls 
and  wall  motion  will  t^e  place  more  easily.  Cyclic  deformation  above  this  critical  temperature, 
will  lead  to  immediate  formation  of  large  cracks,  hence  to  rapid  failure.  Cyclic  deformation 
below  the  critical  temperature  involves  much  smaller  wall  motions,  causing  microcracks  of  a  size 
below  the  detection  limit.  However,  if  a  predeformed  specimen  is  subsequently  deformed  above 
the  critical  temperature,  the  microcracks  will  rapidly  coalesce  to  form  larger  cracks. 

c.  Implications 

★The  irreversibility  which  leads  to  failure  of  PZT  in  cyclic  stressing  is  linked  with  the  motion 
of  ferroelastic  domain  walls.  A  bias  towards  the  direction  of  wall  motion  which  favors  depoling 
results  in  a  ratchet  mechanism.  This  leads  to  grain  boundary  stresses  which  increase  steadily 
until  cracks  form  and  the  material  fails. 

d.  Plans 

★Ferroelastic  domain  wall  motion  will  be  modelled,  below  and  above  the  critical  misfit 
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temperature.  Lifetimes  will  be  calculated  and  linked  to  measurable  parameters.  Critical  questions 
will  be  answered:  specifically,  can  lifetimes  be  lengthened  by  bias  electric  or  stress  field 
stabilization? 

★Confirming  experiments  will  be  suggested  to  White  at  NIST  and  Sadananda  at  NRL. 
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Appendix 


Analytical  Treatment  of  the  Strength  of  Random  Dispersions 


Abstract:  The  methods  of  Boltzmann's  classical  statistical  mechanics  are  employed 
to  calculate  the  yield  strength  in  the  absence  of  thermal  activation  in  the  case 
where  hardening  is  due  to  randomly  dispersed  spherical  particles. 

Introduction 

Ue  seek  to  determine  the  yield  strength  at  0  K  of  a  material  hardened  by 
the  presence  of  spherical  particles  which  act  as  simple  blocks  to  the  motion  of 
dislocations.  The  analysis  is  based  on  a  consideration  of  the  stability  of  a 
sequence  of  lengths  of  a  single  dislocation  which  has  impinged  on  spherical 
particles  in  the  process  of  traversing  a  slip  plane. 

Analysis 

We  assume: 

that  when  acted  on  by  a  uniform  stress  a,  dislocations  are  bent  into  arcs  of 
circles  of  radius  R, 


where  n  is  the  shear  modulus  and  b  the  magnitude  of  the  Burgers  vector; 
that  impenetrable  particles  of  shear  modulus  m  randomly  distributed  and 

provide  restraining  forces,  F,  of  magnitudes  which  are  related  through  the 
breaking  angles  a(a)  to  the  radii  (a)  of  the  circles  of  intersection  between 
particles  and  slip  planes  and  which  are  specified  by 


F  =  a  (a) 


and; 

that  the  yield  stress  is  the  largest  for  which  it  is  possible  to  have  a  stable 
dislocation  join  an  infinite  number  of  particles  along  a  line  which  is 
essentially  straight  and  hence  does  not  cross  itself. 

Towards  the  determination  of  this  stress  we  shall  adapt  the  theoretical 
analysis  first  introduced  by  Hanson  and  Morris  [1]  and  subsequently  modified  and 
employed  by  Labusch  [2]  and  by  Louat  (3)  to  determine  the  yield  strength  due  to 
a  dispersion  of  point  obstacles  all  of  which  provide  the  same  maximum  restraining 
force, 


F  ~  nij^sino-ojijb^  (3) 


This  analysis  rests  on  the  rolling  circle  technique  employed  by  Foreman  and  Makin 
[4]  in  their  numerical  investigation  of  this  hardening.  Thus,  it  was  argued  (see 


fig.  1)  that  a  dislocation  stable  at  a  sequence  of  positions  -  -  ',k-2,  k-1,  will 
be  also  be  stable  at  the  next  position  k  provided  the  particle  associated  with 
the  k+l***  particle  is  centered  anywhere  in  the  region  traversed  by  a  rolling 
circle  as  shown  hatched;  this  region  is  bounded  by  an  angle  0  <  2a.  A  particular 
characteristic  of  such  sequences  of  dislocation  lengths  is  that  they  must  not 
cross  themselves. 

Here,  we  shall  adapt  the  treatment  of  Hanson  and  Morris  [1]  to  the  case 
where  the  strength  of  a  particular  obstacle  depends  on  the  radius  of  the  circle 
of  intersection  between  the  spherical  particle  concerned  and  the  slip  plane. 
To  proceed,  we  note  first  from  fig  (1),  in  the  case  of  point  obstacles,  that, 
following  the  path  of  a  rolling  circle,  the  angle  turned  between  successive 
tangents,  t,^,  is  equal  to  the  angle, which  gives  the  direction  of  the  trailing 
circle  edge  at  the  point  of  rotation  (see  fig  (2)).  It  is  then  readily  seen  that 
in  the  case  of  an  obstacle  with  radius  ,  a,  the  angle  turned  becomes: 

^  =  <j>  (  1  +  a  /  R  )  *2\> 

Then,  if  the  number  of  points  be  N  and  if  the  angle  0  occurs  n;  times  we  have 

(3) 

i-l 

and  the  total  angle  turned  is 

E  $7=  0  (4) 

i-l 

If  N  is  indefinitely  large,  but  ♦  finite,  the  average  value  of  is  zero.  In 
this  case  the  whole  dislocation  line  is  either  straight  or  contains  equal 
infinities  of  complete  clockwise  and  anticlockwise  rotations.  The  latter 
configuration  is  a  situation  of  low  probability  which  we  can  neglect. 

The  procedure  for  the  determination  of  the  thermodynamic  probability  is 
akin  to  that  employed  to  determine  the  Maxwell -Boltzmann  distribution  of  kinetic 
energies  in  an  ideal  monatomic  gas.  Thus,  if  there  are  n  partitions  at  each  of 
N  points  the  thermodynamic  probability  that  the  i*  is  employed  n;  times  is: 


(N  n  flp"' 

p  °  „  '•*  <5> 

n  II  ! 

i-l 

where  p,  is  the  occupational  probability  of  the  i*^  region.  To  determine  p^  we 
divide  the  hatched  area  ,A,  into  n  regions  of  area,  6Aj  (i  =  1  ,n)  each 
characterized  by  a  value  of  say  and  see  that  when  the  concentration  of 
second  phase  particles  is  C  and  provided  the  area  AAj  is  sufficiently  small,  the 
probability  that  a  point  k  lies  in  a  particular  area  AAj  is: 

Here,  we  have  expressed  C  in  terms  of  the  volume  fraction,  f,  of  second  phase. 
Thus,  the  number,  M,  of  particles  per  unit  volume  is  such  that 


2 


-3  «  f  AAj 

=  1-e  -  cAAj  (6) 


Anz^M  _  ^ 
3 


So  that  the  concentration,  C,  of  circles  of  intersection  of  particles  by  a  slip 
plane  is 

C=2r  M  =  . 

2iir* 


On  the  basis  that  the  most  probable  configuration  of  a  system  having  a  large 
number  of  parts  is  that  which  actually  occurs  we  have,  to  determine  this 
configuration,  only  to  maximize  P  subject  to  the  constraints  offered  by  (3) 
and(4).  To  proceed  we  take  advantage  of  the  facts  that  N  and  n^  are  both  large 
to  employ  Stirling's  approximation  for  the  factorial  and  then  set  the  variation 
in  In  P  to  zero.  We  find: 


while  (3)  and(4)  become: 


(7) 

i-i 


E  =  0 .  (8) 

Whence  introducing  the  Lagrangian  multipliers  X  and  y  we  have 

In  Pj-  In  n_j+  Y  +  0 


so  that 


N  pj©*^ 


n 

£Pi  e 

i-l 


»4i 


(9) 


Then  from  (4)  and  (8)  we  have  so  that  we  find  on  substituting  for  O;  in  the 
logarithm  of  P  and  after  some  manipulation: 


3 


NL  pj  =  *  5, 


where 

=  5.  (10) 


P  =  (11) 

P  *  s'*  e"'*.  (11) 

Then  invoking  (6)  and  (9)  and  proceeding  to  the  limit  of  small  quantities  S  A 
,  (10)  becomes: 


P  =  [  cj^e*-^  dA]*^ 


We  have  now  to  determine  S  and  X  .  Here,  where  A  is  the  area  indicated  in 
fig.(l)  we  have,  for  1/2R«1: 

2It  r(2*-aic8ini  /2R) 

S  =  cj  (1  +r  )  dl  J 

0  T(arc8in  i  /2R) 


^  Jg2r  (g28Ar-i)  (n-e-n^r) 

X  (XD2+1 


(12) 


The  determination  of  X  is  assisted  by  the  realization  that 


d  Ins  =  1  ^  ^0 

dX  ~  s  dX  ~  f:'  s  ~  n" 


as  iV  -*  «> 


(13) 


and  X  must  be  such  that 


dS 

dX 


=  0. 


This  question  can  only  be  resolved  numerically.  In  the  special  and  unrealistic 
case  in  which  a  «  1  we  find,  independent  of  the  value  of  T,  that  X  *  1.41  /  a. 
In  comparison  with  the  case  treated  previously  [1]  there  are  additional 
complications:  the  obstacles  have  finite  radius  and;  because  the  radii  of 
individual  obstacles  differ  one  from  the  another,  they  have  variable  strengths. 
We  shall  now  explore  the  nature  of  this  variation.  Bacon  ,Kocks  and  Scattergood 
[4]  have  shown  through  numerical  studies  that  the  Orowan  stress, r,  for  a  linear 
sequence  of  impenetrable  circular  blocks  to  dislocation  motion  having  radius  a 
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and  Inter-center  spacing  L  is  given  by: 


T 


= 

4ll 


In(^) 

_ D 

(L-2a) 


where, 


1  _  1 
X  ~  2a 


(14) 


Thus,  the  specific  strength  of  these  obstacles  is  variable  and  denoted  by  the 
term: 


ln(x) 

4n 


=  sin  a. 


(15) 


Here,  the  obstacles  are  not  distributed  along  a  straight  line  nor  are  they 
regularly  spaced.  Additionally  their  strengths  are  variable.  The  consequences 
of  the  first  two  particulars  have  been  considered  in  the  foregoing,  we  now 
examine  the  consequences  of  the  variation  In  specific  strength. 

To  progress  we  suppose  that  the  particles  are  spherical  of  radius  r  and 
that  the  concentration  of  the  particulate  phase  is  f.  We  then  find  that  the  mean 
spacing  between  particle  centers  In  a  slip  plane  Is 


Again,  the  mean  value  of  the  radius  of  Intersection  of  such  a  plane  with 
the  particles  Is  readily  found  to  be: 


Then  towards  the  evaluation  of  the  quantity  defined  in  (15)  we  follow 
Friedel  and  suppose  that  at  yield,  dislocation  are  bent  into  arcs  of  circles  of 
radius  R  such  that  the  areas  swept  by  the  dislocation  In  circular  segments 
between  obstacles  contain  about  one  obstacle.  Taking  the  critical  angle  turned 
by  the  dislocation  at  an  obstacle  as  29,  the  area  of  a  segment  Is: 


4j;^e^ 

3 


Thence  we  see  that 
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^  (L-2a)^e.f  _  2icr^ 

3  3  ”  3 


so  that 


I.-2a 

a 


and 


L  =  a  (2+-, 


N  JfB 


Thus,  the  mean  value  of  the  quantity  X/b  defined  in  (14)  is 

2r 

b 

and  provided,  as  is  usually  the  case,  f  ^  »  0.1, 

(15)  becomes 


4n 


sina 


Importantly  in  the  circumstances  considered  here,  the  strength  of  a  particular 
obstacle  of  radius  a  is  denoted  by 


sin 


4n 


(16) 


a<~60<>. 


An  immediate  consequence  of  these  conclusions  is  that  the  radius  R  »a  since  it 
is  significantly  larger  than  L  unless  0  >  60*. 

We  are  now  in  a  position  to  consider  the  consequences  of  the  simultaneous 
existence  of  obstacles  of  differing  strengths.  Reference  to  (12)  and  the 
realization  that  the  analysis  which  proceeds  it,  is  essentially  linear  in  the 
concentration  of  obstacles,  shows  that  we  may  define  a  parameter  Sj  for  each 
group  of  obstacles  characterised  by  a  strength  otj  and  a  concentration  Cj  and  then 
consider  the  quantity: 

which  in  the  limit  of  a  continuous  distribution  of  obstacle  sizes  becomes: 
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R^JlCAe^‘^^-1) 

_ _ 

^(X^T+l) 


— >  fg(tt)  (e^‘^^-l)da.  (17) 

X (k^T+1)  { 

where  is  the  maximum  strength  of  any  obstacle,  and  Z„  is  the  concentration 
of  obstacles  of  strength  at. 

In  the  case  under  consideration  where  the  particles  are  spherical  and  all 
of  the  same  radius,  r,  it  is  easily  shown  that  the  fraction,  f(a)  of  circles  of 
intersection  of  having  radii  in  the  range  a  to  a  +5a  is 

C{a)ba  =  - —  ha. 

r  ^r^-a^ 


and  substituting  for  a  in  terms  of  a,  (17)  becomes: 


So  =  fr- 


■/ 


ll*a/R)X 

(23)  21.  _i 

D 


4firX)  J  a^(i+a/R)*l) 


da, 


where  2fr  is  the  concentration  of  obstacles  in  a  plane. 

To  progress,  we  observe  that  r/R  is  small  and  that  the  integral,  I,  lies 
between  lower  and  upper  bounds.  Thus, 


fR‘ 


/ 


2ax  ali 


da 


<5„ 


fR2(^)  2mR  ^  (^)  2* 

<— ^ - J-r^da. 


8X{X^+1) 


(18) 


The  integral  I  can  be  expressed  in  terms  of  the  Beta  function.  Thus,  we  find: 
The  interpretation  of  this  result  is  assisted  by  the  use  of  Stirling's 
approximation: 

Whence  we  find  that  (19)  is: 
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(19) 


c 

2VA,+2it  X(X2+1) 
<^o 

i;2  c  r 

2vT+2W  X(A,2+1) 


Inspection  shows  that  these  bounds  are  nearly  coincident  so  that  we  may 
concentrate  our  attention  on  the  simpler  of  the  two. 

Towards  the  determination  of  \  we  have  as  indicated  earlier: 

d  ln5o  ^  X  ^  -  S 
dX  Sq  dX  '  So  N' 

a  quantity  which  tends  to  0  as  N  tends  to  infinity. 
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Accordingly  differentiating  with  regard  to  X  we  find  that  we  require 

In^  -  -1=0. 

2ii  b  X+2ii  X^+l  A, 

In  the  range  of  interest  in  which 

102  i  ^  i  10® 


we  find  that 

22X^1. 


Finally  we  require  as  a  condition  for  yield  that  So  =  1.  Accordingly 
setting 


in  (12)  and  substituting  for  C  in  terms  of  f  we  have 

— 

3^2Jb2/(^)  2. 

^ - 

471  r Vl‘^2n  X(X^+1) 

Since  X  is  largest  when  r  is  least  we  find  that  there  is  little  variation 
in  the  term  (2r/b)'’'^"  with  variation  in  X.  Overall,  for  a  constant  value  of  f, 
the  yield  stress  varies  nearly  with  the  inverse  of  the  particle  size.  It  may  be 
noted  that  Freidel's  [4]  conclusion  that  stress  varies  with  the  square  root  of 
concentration  is  retained. 

Conclusions: 

An  analytical  analysis  which  proceeds  on  the  same  lines  as  those  previously 
employed  to  evaluate  the  hardening  effect  of  a  dispersion  of  point  obstacles  has 
been  used  to  deal  with  the  case  where  the  particles  have  finite  radius. 
Qualitatively,  the  principle  conclusions  are  that  the  predicted  hardening  should 
increase  with  the  square  root  of  the  concentration  of  particulate  phase  and 
approximately  with  inverse  of  the  particle  diameter. 
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